Comparative analysis of Wolbachia maternal transmission and localization in host ovaries.

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Coevolution with toxic prey produces functional trade-offs in sodium channels of predatory snakes.

Seemingly unrelated traits often share the same underlying molecular mechanisms, potentially generating a pleiotropic relationship whereby selection shaping one trait can simultaneously compromise another. While such functional trade-offs are expected to influence evolutionary outcomes, their actual relevance in nature is masked by obscure links between genotype, phenotype, and fitness. Here, we describe functional trade-offs that likely govern a key adaptation and coevolutionary dynamics in a predator-prey system. Several garter snake (Thamnophis spp.) populations have evolved resistance to tetrodotoxin (TTX), a potent chemical defense in their prey, toxic newts (Taricha spp.). Snakes achieve TTX resistance through mutations occurring at toxin-binding sites in the pore of snake skeletal muscle voltage-gated sodium channels (NaV1.4). We hypothesized that these mutations impair basic NaV functions, producing molecular trade-offs that should ultimately scale up to compromised organismal performance. We investigate biophysical costs in two snake species with unique and independently evolved mutations that confer TTX resistance. We show electrophysiological evidence that skeletal muscle sodium channels encoded by toxin-resistant alleles are functionally compromised. Furthermore, skeletal muscles from snakes with resistance genotypes exhibit reduced mechanical performance. Lastly, modeling the molecular stability of these sodium channel variants partially explains the electrophysiological and muscle impairments. Ultimately, adaptive genetic changes favoring toxin resistance appear to negatively impact sodium channel function, skeletal muscle strength, and organismal performance. These functional trade-offs at the cellular and organ levels appear to underpin locomotor deficits observed in resistant snakes and may explain variation in the population-level success of toxin-resistant alleles across the landscape, ultimately shaping the trajectory of snake-newt coevolution.

Distinct Wolbachia localization patterns in oocytes of diverse host species reveal multiple strategies of maternal transmission.

A broad array of endosymbionts radiate through host populations via vertical transmission, yet much remains unknown concerning the cellular basis, diversity, and routes underlying this transmission strategy. Here, we address these issues, by examining the cellular distributions of Wolbachia strains that diverged up to 50 million years ago in the oocytes of 18 divergent Drosophila species. This analysis revealed 3 Wolbachia distribution patterns: (1) a tight clustering at the posterior pole plasm (the site of germline formation); (2) a concentration at the posterior pole plasm, but with a significant bacteria population distributed throughout the oocyte; and (3) a distribution throughout the oocyte, with none or very few located at the posterior pole plasm. Examination of this latter class indicates Wolbachia accesses the posterior pole plasm during the interval between late oogenesis and the blastoderm formation. We also find that 1 Wolbachia strain in this class concentrates in the posterior somatic follicle cells that encompass the pole plasm of the developing oocyte. In contrast, strains in which Wolbachia concentrate at the posterior pole plasm generally exhibit no or few Wolbachia in the follicle cells associated with the pole plasm. Taken together, these studies suggest that for some Drosophila species, Wolbachia invade the germline from neighboring somatic follicle cells. Phylogenomic analysis indicates that closely related Wolbachia strains tend to exhibit similar patterns of posterior localization, suggesting that specific localization strategies are a function of Wolbachia-associated factors. Previous studies revealed that endosymbionts rely on 1 of 2 distinct routes of vertical transmission: continuous maintenance in the germline (germline-to-germline) or a more circuitous route via the soma (germline-to-soma-to-germline). Here, we provide compelling evidence that Wolbachia strains infecting Drosophila species maintain the diverse arrays of cellular mechanisms necessary for both of these distinct transmission routes. This characteristic may account for its ability to infect and spread globally through a vast range of host insect species.

Conspicuous coloration of toxin-resistant predators implicates additional trophic interactions in a predator-prey arms race.

Antagonistic coevolution between natural enemies can produce highly exaggerated traits, such as prey toxins and predator resistance. This reciprocal process of adaptation and counter-adaptation may also open doors to other evolutionary novelties not directly involved in the phenotypic interface of coevolution. We tested the hypothesis that predator-prey coevolution coincided with the evolution of conspicuous coloration on resistant predators that retain prey toxins. In western North America, common garter snakes (Thamnophis sirtalis) have evolved extreme resistance to tetrodotoxin (TTX) in the coevolutionary arms race with their deadly prey, Pacific newts (Taricha spp.). TTX-resistant snakes can retain large amounts of ingested TTX, which could serve as a deterrent against the snakes' own predators if TTX toxicity and resistance are coupled with a conspicuous warning signal. We evaluated whether arms race escalation covaries with bright red coloration in snake populations across the geographic mosaic of coevolution. Snake colour variation departs from the neutral expectations of population genetic structure and covaries with escalating clines of newt TTX and snake resistance at two coevolutionary hotspots. In the Pacific Northwest, bright red coloration fits an expected pattern of an aposematic warning to avian predators: TTX-resistant snakes that consume highly toxic newts also have relatively large, reddish-orange dorsal blotches. Snake coloration also seems to have evolved with the arms race in California, but overall patterns are less intuitively consistent with aposematism. These results suggest that interactions with additional trophic levels can generate novel traits as a cascading consequence of arms race coevolution across the geographic mosaic.

Temperature effects on cellular host-microbe interactions explain continent-wide endosymbiont prevalence.

Endosymbioses influence host physiology, reproduction, and fitness, but these relationships require efficient microbe transmission between host generations to persist. Maternally transmitted Wolbachia are the most common known endosymbionts,1 but their frequencies vary widely within and among host populations for unknown reasons.2,3 Here, we integrate genomic, cellular, and phenotypic analyses with mathematical models to provide an unexpectedly simple explanation for global wMel Wolbachia prevalence in Drosophila melanogaster. Cooling temperatures decrease wMel cellular abundance at a key stage of host oogenesis, producing temperature-dependent variation in maternal transmission that plausibly explains latitudinal clines of wMel frequencies on multiple continents. wMel sampled from a temperate climate targets the germline more efficiently in the cold than a recently differentiated tropical variant (∼2,200 years ago), indicative of rapid wMel adaptation to climate. Genomic analyses identify a very narrow list of wMel alleles-most notably, a derived stop codon in the major Wolbachia surface protein WspB-that underlie thermal sensitivity of cellular Wolbachia abundance and covary with temperature globally. Decoupling temperate wMel and host genomes further reduces transmission in the cold, a pattern that is characteristic of host-microbe co-adaptation to a temperate climate. Complex interactions among Wolbachia, hosts, and the environment (GxGxE) mediate wMel cellular abundance and maternal transmission, implicating temperature as a key determinant of Wolbachia spread and equilibrium frequencies, in conjunction with Wolbachia effects on host fitness and reproduction.4,5 Our results motivate the strategic use of locally selected wMel variants for Wolbachia-based biocontrol efforts, which protect millions of individuals from arboviruses that cause human disease.6.

Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts.

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Nav) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Nav gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Nav paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Nav gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Nav1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Nav1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

Pervasive effects of Wolbachia on host activity.

Heritable symbionts have diverse effects on the physiology, reproduction and fitness of their hosts. Maternally transmitted Wolbachia are one of the most common endosymbionts in nature, infecting about half of all insect species. We test the hypothesis that Wolbachia alter host behaviour by assessing the effects of 14 different Wolbachia strains on the locomotor activity of nine Drosophila host species. We find that Wolbachia alter the activity of six different host genotypes, including all hosts in our assay infected with wRi-like Wolbachia strains (wRi, wSuz and wAur), which have rapidly spread among Drosophila species in about the last 14 000 years. While Wolbachia effects on host activity were common, the direction of these effects varied unpredictably and sometimes depended on host sex. We hypothesize that the prominent effects of wRi-like Wolbachia may be explained by patterns of Wolbachia titre and localization within host somatic tissues, particularly in the central nervous system. Our findings support the view that Wolbachia have wide-ranging effects on host behaviour. The fitness consequences of these behavioural modifications are important for understanding the evolution of host-symbiont interactions, including how Wolbachia spread within host populations.

Pervasive Effects of Wolbachia on Host Temperature Preference.

Heritable symbionts can modify a range of ecologically important host traits, including behavior. About half of all insect species are infected with maternally transmitted Wolbachia, a bacterial endosymbiont known to alter host reproduction, nutrient acquisition, and virus susceptibility. Here, we broadly test the hypothesis that Wolbachia modifies host behavior by assessing the effects of eight different Wolbachia strains on the temperature preference of six Drosophila melanogaster subgroup species. Four of the seven host genotypes infected with A-group Wolbachia strains (wRi in Drosophila simulans, wHa in D. simulans, wSh in Drosophila sechellia, and wTei in Drosophila teissieri) prefer significantly cooler temperatures relative to uninfected genotypes. Contrastingly, when infected with divergent B-group wMau, Drosophila mauritiana prefers a warmer temperature. For most strains, changes to host temperature preference do not alter Wolbachia titer. However, males infected with wSh and wTei tend to experience an increase in titer when shifted to a cooler temperature for 24 h, suggesting that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results indicate that Wolbachia modifications to host temperature preference are likely widespread, which has important implications for insect thermoregulation and physiology. Understanding the fitness consequences of these Wolbachia effects is crucial for predicting evolutionary outcomes of host-symbiont interactions, including how Wolbachia spreads to become common.IMPORTANCE Microbes infect a diversity of species, influencing the performance and fitness of their hosts. Maternally transmitted Wolbachia bacteria infect most insects and other arthropods, making these bacteria some of the most common endosymbionts in nature. Despite their global prevalence, it remains mostly unknown how Wolbachia influence host physiology and behavior to proliferate. We demonstrate pervasive effects of Wolbachia on Drosophila temperature preference. Most hosts infected with A-group Wolbachia prefer cooler temperatures, whereas the one host species infected with divergent B-group Wolbachia prefers warmer temperatures, relative to uninfected genotypes. Changes to host temperature preference generally do not alter Wolbachia abundance in host tissues, but for some A-group strains, adult males have increased Wolbachia titer when shifted to a cooler temperature. This suggests that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results help elucidate the impact of endosymbionts on their hosts amid the global Wolbachia pandemic.

The geographic mosaic of arms race coevolution is closely matched to prey population structure.

Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of "trait remixing" in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

Environmental and Genetic Contributions to Imperfect wMel-Like Wolbachia Transmission and Frequency Variation.

Maternally transmitted Wolbachia bacteria infect about half of all insect species. They usually show imperfect maternal transmission and often produce cytoplasmic incompatibility (CI). Irrespective of CI, Wolbachia frequencies tend to increase when rare only if they benefit host fitness. Several Wolbachia, including wMel that infects Drosophila melanogaster, cause weak or no CI and persist at intermediate frequencies. On the island of São Tomé off West Africa, the frequencies of wMel-like Wolbachia infecting Drosophila yakuba (wYak) and Drosophila santomea (wSan) fluctuate, and the contributions of imperfect maternal transmission, fitness effects, and CI to these fluctuations are unknown. We demonstrate spatial variation in wYak frequency and transmission on São Tomé. Concurrent field estimates of imperfect maternal transmission do not predict spatial variation in wYak frequencies, which are highest at high altitudes where maternal transmission is the most imperfect. Genomic and genetic analyses provide little support for D. yakuba effects on wYak transmission. Instead, rearing at cool temperatures reduces wYak titer and increases imperfect transmission to levels observed on São Tomé. Using mathematical models of Wolbachia frequency dynamics and equilibria, we infer that temporally variable imperfect transmission or spatially variable effects on host fitness and reproduction are required to explain wYak frequencies. In contrast, spatially stable wSan frequencies are plausibly explained by imperfect transmission, modest fitness effects, and weak CI. Our results provide insight into causes of wMel-like frequency variation in divergent hosts. Understanding this variation is crucial to explain Wolbachia spread and to improve wMel biocontrol of human disease in transinfected mosquito systems.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis).

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Nav proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Nav1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

Large-effect mutations generate trade-off between predatory and locomotor ability during arms race coevolution with deadly prey.

Adaptive evolution in response to one selective challenge may disrupt other important aspects of performance. Such evolutionary trade-offs are predicted to arise in the process of local adaptation, but it is unclear if these phenotypic compromises result from the antagonistic effects of simple amino acid substitutions. We tested for trade-offs associated with beneficial mutations that confer tetrodotoxin (TTX) resistance in the voltage-gated sodium channel (NaV1.4) in skeletal muscle of the common garter snake (Thamnophis sirtalis). Separate lineages in California and the Pacific Northwest independently evolved TTX-resistant changes to the pore of NaV1.4 as a result of arms race coevolution with toxic prey, newts of the genus Taricha. Snakes from the California lineage that were homozygous for an allele known to confer large increases in toxin resistance (NaV1.4LVNV) had significantly reduced crawl speed compared to individuals with the ancestral TTX-sensitive channel. Heterologous expression of native snake NaV1.4 proteins demonstrated that the same NaV1.4LVNV allele confers a dramatic increase in TTX resistance and a correlated decrease in overall channel excitability. Our results suggest the same mutations that accumulate during arms race coevolution and beneficially interfere with toxin-binding also cause changes in electrophysiological function of the channel that may affect organismal performance. This trade-off was only evident in the predator lineage where coevolution has led to the most extreme resistance phenotype, determined by four critical amino acid substitutions. If these biophysical changes also translate to a fitness cost-for example, through the inability of T. sirtalis to quickly escape predators-then pleiotropy at this single locus could contribute to observed variation in levels of TTX resistance across the mosaic landscape of coevolution.

Convergent adaptation to dangerous prey proceeds through the same first-step mutation in the garter snake Thamnophis sirtalis.

Convergent phenotypes often result from similar underlying genetics, but recent work suggests convergence may also occur in the historical order of substitutions en route to an adaptive outcome. We characterized convergence in the mutational steps to two independent outcomes of tetrodotoxin (TTX) resistance in separate geographic lineages of the common garter snake (Thamnophis sirtalis) that coevolved with toxic newts. Resistance is largely conferred by amino acid changes in the skeletal muscle sodium channel (NaV 1.4) that interfere with TTX-binding. We sampled variation in NaV 1.4 throughout western North America and found clear evidence that TTX-resistant changes in both lineages began with the same isoleucine-valine mutation (I1561V) within the outer pore of NaV 1.4. Other point mutations in the pore, shown to confer much greater resistance, accumulate later in the evolutionary progression and always occur together with the initial I1561V change. A gene tree of NaV 1.4 suggests the I1561V mutations in each lineage are not identical-by-decent, but rather they arose independently. Convergence in the evolution of channel resistance is likely the result of shared biases in the two lineages of T. sirtalis-only a few mutational routes can confer TTX resistance while maintaining the conserved function of voltage-gated sodium channels.

Toxicity and population structure of the Rough-Skinned Newt (Taricha granulosa) outside the range of an arms race with resistant predators.

Species interactions, and their fitness consequences, vary across the geographic range of a coevolutionary relationship. This spatial heterogeneity in reciprocal selection is predicted to generate a geographic mosaic of local adaptation, wherein coevolutionary traits are phenotypically variable from one location to the next. Under this framework, allopatric populations should lack variation in coevolutionary traits due to the absence of reciprocal selection. We examine phenotypic variation in tetrodotoxin (TTX) toxicity of the Rough-Skinned Newt (Taricha granulosa) in regions of allopatry with its TTX-resistant predator, the Common Garter Snake (Thamnophis sirtalis). In sympatry, geographic patterns of phenotypic exaggeration in toxicity and toxin-resistance are closely correlated in prey and predator, implying that reciprocal selection drives phenotypic variation in coevolutionary traits. Therefore, in allopatry with TTX-resistant predators, we expect to find uniformly low levels of newt toxicity. We characterized TTX toxicity in northwestern North America, including the Alaskan panhandle where Ta. granulosa occur in allopatry with Th. sirtalis. First, we used microsatellite markers to estimate population genetic structure and determine if any phenotypic variation in toxicity might be explained by historical divergence. We found northern populations of Ta. granulosa generally lacked population structure in a pattern consistent with northern range expansion after the Pleistocene. Next, we chose a cluster of sites in Alaska, which uniformly lacked genetic divergence, to test for phenotypic divergence in toxicity. As predicted, overall levels of newt toxicity were low; however, we also detected unexpected among- and within-population variation in toxicity. Most notably, a small number of individuals contained large doses of TTX that rival means of toxic populations in sympatry with Th. sirtalis. Phenotypic variation in toxicity, despite limited neutral genetic divergence, suggests that factors other than reciprocal selection with Th. sirtalis likely contribute to geographic patterns of toxicity in Ta. granulosa.